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We contrast the evolution of gravitational binary mergers for the two cases of a black 

m ; 

C^) , hole and a normal neutron star, and a black hole and a self-bound strange quark matter 

O 

. star. In both cases inspiral continues until the Roche limit is reached, at which point it 

£>y is expected that stable mass transfer to the balck hole ensues. Whereas a neutron star 

would then outspiral, the strange quark matter star barely moves out at all. Eventually, 



the strange quark star loses all its mass to the black hole, as compared to the neutron star 
whose outspiral continues until stable mass transfer terminates and inspiral resumes. These 
scenarios result in distinctly different gravitational wave signatures. 
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r-j '. Several proposals concerning the physical state and the internal constitution of matter at supra- 

Or 



nuclear densities have been put forth (see |lj,l2j ^ or recen t accounts). Figure Q shows many exciting 

O ; 

+3 , possibilities for the composition of compact stars including (1) strangeness-bearing matter in the 

form of hyperons, kaons, or quarks, (2) Bose (pion or kaon) condensed matter, and (3) so-called 



self-bound strange quark matter (SQM). Fermions, whether they are in the form of baryons or 



deconfined quarks, are expected to additionally exhibit superfluidity and/or superconductivity. In 
the case of quarks, this could occur either as two-flavor superconductivity (2SC) or in the Color- 
Flavor-Locked (CFL) phase in which all three flavors (up, down, and strange) participate. 

The possibility of such exotic phases brings attendant changes to the predictions of maximum 
masses and radii, since the presence of multiple components or new phases of matter generally 
lessens the pressure for a given energy density. The fact that compact objects are the only objects 
in which such phases could occur underscores the importance of precise determinations of basic 
observables: 

(i) masses and radii [J]. Simply put, a precise determination of the mass and radius of the same 
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FIG. 1: A schematic diagram of the possible phases in a compact star (Courtesy: A.W. Steiner). 

neutron star would not only be a first achievement for observational astronomy, but would also 
hold the promise of delineating Quantum Chromo-Dynamics (QCD) in regimes of high baryon 
density heretofore not possible. 

(ii) surface temperatures versus age [J. Detections of photons from cooling neutron stars could 
constrain the pairing gaps of baryons and the star's mass. However, establishing the sizes of quark 
gaps or the presence of Goldstone bosons in the CFL phase of quarks in a hybrid star will be 
difficult ^, I4 1 because the layers of normal matter surrounding the quark core continues to cool 



through significantly more rapid processes. In contrast, the mean energy of emitted photons from 



lan that from a normal 
. Due to its distinctive 



the bare surface of a strange quark star would be significantly larger t 
cooling neutron star (30 < E/keV < 500 versus 0.1 < E/keV < 2.5) j|i 
spectrum and time evolution, such an observation would constitute an unmistakable detection of 
a strange quark star and shed light on color superconductivity at "stellar" densities, 
(hi) anomalous behavior in the spin-down rates of ms radio pulsars [f| . Possible hysteresis (reversal) 
in the normal spindown rate of pulsars would signal the appearance of a mixed-phase of matter 
containing normal hadronic matter with more compressible Bose condensed or quark matter, 
(iv) neutrino luminosities from future galactic core collapse supernovae. The main new feature that 
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has emerged from studies of neutron stars at birth is the possible metastability and subsequent 
collapse to a black hole of a proto-neutron star containing quark matter, or other types of matter 
including hyperons or a Bose condensate, which could be observable in the neutrino signal 0. 

Examples of ongoing and planned observations of solar mass compact objects that could shed 
light on QCD at high baryon density include 

(i) multi-wavelength photon observations with the HST, Chandra, XMM, Integral, etc., 

(ii) spectral and temporal studies of supernova neutrinos with the SuperK, SNO, UNO, etc., and 

(iii) gravity wave detections from coalescing compact object binaries though LIGO, VIRGO, etc. 
In this work, we contrast the evolution of binary star mergers for two distinct cases: 

(1) A black hole (BH) and a normal star. For the discussion at hand, this refers to a star with a 
surface of normal matter in which the pressure vanishes at vanishing baryon density. The interior 
of the star, however, may contain any or a combination of the many exotica [3] depicted in Fig. ^ 

(2) A BH and a self-bound star which is exemplified by Witten's Strange-Quark Matter (SQM) 
star 0|; see for a review. Such a star has a bare quark matter surface in which the pressure 
vanishes at a finite but supra- nuclear baryon density [2J|. 

Prototypes of the mass versus central baryon density (n c ) and mass versus radius for these 
cases are shown in Fig. |5J Quantitative variations from these generic behaviors can be caused by 
uncertainties in the underlying strong interaction models (see the compendium of results in Fig. 2 
of Ref. Qualitative differences in the outcomes of mergers with a black hole emerge, however, 
because of the gross differences in the mass-radius diagram. 
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FIG. 2: (a) Mass versus central density, and (b) mass versus radius, for normal and self-bound stars. 
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A normal star and a self-bound star represent two quite different possibilities (see Fig. |3J) for 
the quantity 



dlnR I < for a normal neutron star (NS) 
dlnM I > o for a se if. bounc i SQM star 



where M and R are the star's mass and radius, respectively. For small to moderate mass self-bound 
stars, R oc M 1 / 3 so that a = 1/3; only for configurations approaching the maximum mass does a 
turn negative. 
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FIG. 3: Comparisons of a and radii R for the normal star (NS) and self-bound strange quark matter (SQM) 
star as functions of mass. 

Our objective is to explore the astrophysical consequences of these distinctive behaviors in R 
versus M as they affect mergers with a black hole. Note that a is intimately connected with 
the dense matter equation of state (EOS), since there exists a one-to-one correspondence between 
R(M) and P(ns), where P is the pressure and hb is the baryon density. Gravitational mergers 
in which a compact star loses its mass (either to a companion star or to an accretion disk) during 
evolution is one of the rare examples in which the R versus M (or equivalently, P versus rig) 
relationship of the same star is sampled. Although we focus here on the coalescence of a compact 
star with a BH, the theoretical formalism and our principal findings apply also to mergers in which 
both objects are compact stars. 
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II. THE MERGER OF A COMPACT STAR WITH A LOW-MASS BLACK HOLE 

The general problem of the origin and evolution of systems containing a neutron star and a 
black hole was first detailed by Lattimer & Schramm [12I (lffi . Compact binaries form naturally as 
the result of evolution of massive stellar binaries (e.g., The estimated lower mass limit for 

supernovae (which produce neutron stars or black holes) is approximately 8 Mq. Observationally, 
the number of binaries formed within a given logarithmic separation is approximately constant, 
but the relative mass distributions are uncertain. There is some indication that the distribution 
in binary mass ratios might also be flat. Most progenitor systems do not survive the first, more 
massive star becoming a supernova. In the absence of a "kick" velocity from the explosion, it is 
easily found that the loss of more than half of the mass from the system will unbind it. However, 
the fact that pulsars are observed to have mean velocites in excess of a few hundred km/s implies 
that neutron stars are usually produced with large kicks. In the cases that the kick, which is 
thought to be randomly directed, opposes the star's orbital velocity, the chances that the post- 
supernova binary remains intact increases. In addition, the orbital separation in a surviving binary 
will be reduced significantly. Subsequent evolution then progresses to the supernova explosion of 
the companion. Those systems that survive the second explosion should both have greatly reduced 
separations and orbits with high eccentricity. 

Gravitational radiation causes the binary's orbit to decay such that a system with masses 
Mi and M2 with initial semimajor axes a satisfying 

a < 2.8[M 1 M 2 (M 1 + M 2 )/M|] 1 / 4 R , (2) 

will fully decay within the age of the Universe (~ 10 10 yr). This limit is for circular orbits; those 
with highly eccentric orbits will decay much faster Q- Ref. Q argued that (1) mergers of neutron 
stars and black holes coupled with the subsequent ejection of a few percent of the neutron star's 
mass, could easily account for all the r-process nuclei in the cosmos, and (2) compact object binary 
mergers could be associated with gamma-ray bursts (see also Jjfl ). 

When the less massive inspiralling compact star reaches its Roche limit (see Fig. H}, mass can 
be stripped from it. The radius of the compact object will quickly adjust to its new mass. If 
the radius increases more quickly than the Roche limiting radius, mass transfer to the BH will be 
stable, and the inspiral will be halted due to angular momentum conservation. The classical Roche 
limit is based upon an incompressible fluid of density p and mass Mi in orbit about a mass M±. 





M 1 = 1.5M Sun , M 2 = 3.5M Sun , 


a = 25 M Sun 
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FIG. 4: Equipotential contours of Newtonian gravitational plus centrifugal potential of a compact star and 
a black hole binary. The inner dashed curve shows the Roche lobe (Courtesy: S. Ratkovic). 

In Newtonian gravity, this limit is 

RRochcNewt = (A/1/0.09017TP) 1 / 3 = 19.2(Mi /M 0/ oi 5 ) 1/3 km , (3) 

where pis = p/10 15 g cm~ 3 . Using general relativity, Fishbone found that the number 
0.0901 in Eq. Q becomes 0.0664, even for rotating BHs. In geometrized units, R-Roche/Mi = 
13(14. 4)(M^pi5/Mq) -1 / 3 , where the numerical coefficient refers to the Newtonian (GR) case. In 
other words, if the neutron star's mean density is p\§ = 1, the Roche limit is encountered beyond 
the last stable orbit (R = 6Mi for a non-rotating BH) if Mi < 5.9 M , leading to mass overflow 
and mass transfer. And, as now discussed, the mass transfer may proceed stably for some consid- 
erable time. This would lengthen the lifetime the BH would accrete matter from its companion, 
which, in itself, could be an observational signature. 

III. EVOLUTION OF MERGERS 

The final evolution of a compact binary is now discussed (see Fig. for a schematic illustration). 
Define q = 7ti cs /Mbh, A 4 = m cs^BH /M ", and M = Mbh + mcs, where m cs and Mbh are the 
compact star (normal star or self-bound SQM star) and black hole masses, respectively. The 



orbital angular momentum is 



J 2 = Gf/Ma = GM 3 aq 2 /(l + qf . 
We can employ Paczynski's jlll ] formula for the Roche radius of the secondary: 

R e /a = 0A6[q/(l + q)} 1/3 , 

or a better fit by Eggleton j^ : 

Ri/a = 0.49[0.6 + q~ 2/3 ln(l + g 1/3 )] _1 . 



(4) 



(5) 



(6) 



The orbital separation a during mass transfer is obtained by setting Ri = R, the compact star 
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FIG. 5: A schematic illustration of a compact star merger with a black hole (Courtesy: M.W. Carmcll). 



radius. For stable mass transfer, the star's radius has to increase more quickly than the Roche 
radius as mass is transferred 



251 ] . Thus, we must have, using Paczynski's formula, 



dlnR din Re din a 1 

= a > = 1 — 

dlnm cs d\nm cs dlnm cs 3 

for stable mass transfer. If the mass transfer is conservative j^fj], j = Jew, where 

32 G 1 ! 2 fi 2 M 5 / 2 32 G 7 / 2 q 2 M 9 / 2 



Jew 



and 



5 c 5 a 7 / 2 5 c 5 {l+qfa 7 / 2 



J a q(l 



J 2a <?(! + 



(7) 



(8) 



(9) 
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where the dots on various symbols denote time derivatives. This leads to 




32 G 3 q 2 M 3 
1 7(l + ?)a 4 ' 



(10) 



Since m cs < Mbh, q < 0, and the condition for stable mass transfer is simply 



q < 5/6 + a/2. 



(11) 



This condition is achieved in all normal neutron star - BH or self-bound star - BH binaries. In 
both cases, the luminosities of gravity wave emission can be estimated using 



In addition to Eq. (1111) . a number of other conditions must hold for stable mass transfer to 

n n 

occur (c/. |2J|, and references therein). The discussion in |21|] is relevant to the stellar merger of a 
self-bound star as well. We now compare the mergers of normal and self-bound stars with a BH. For 
definiteness, we choose m ns = msQM = 1-5M and Mbh = 3.5M Q . Results for other choices are 
straightforward to obtain. Fig. (left panel) shows the time development of the orbital separation 
a and the compact star's mass and radius during the stable mass transfer (outspiral) phases. The 
time evolutions during stable mass transfer are obtained from Eq. (|10ft . using m cs = qM/(l + q) 2 . 

In the case of SQM, the inspiral continues to smaller separations since Rsqm < Rns- As a 
consequence, mass transfer is accelerated and the duration of the stable mass transfer outspiral 
phase is shortened considerably. In the SQM case the star barely spirals out since 1/3 — a ~ (see 
Eq. ©). The star eventually loses all its mass. In contrast, for a neutron star 1/3 — a > and the 
star outspirals significantly, until a ~ —5/3 when stable mass transfer terminates (Eq. 1)11(1 ). This 
occurs above the minimum stable neutron star mass, ~ 0.08 M , so the remnant mini-neutron 
star will resume spiralling in, concluding in a final merger. In either case, if stable mass transfer 
occurs, it proceeds for much longer than an orbital period, perhaps up to a few tenths of a second. 
The merger of neutron star and self-bound star with a BH will significantly differ in their gravity 
wave signatures (shown in the right panel of Fig. EJ). The self-bound star case will have a shorter 
timescale of emission, with probably higher peak luminosity, and will lack a final increase due to 
resumption of inspiral in the NS case. 

An additional contrast between mergers involving only NSs or only SQM stars (but not BHs) 
is that the condition for stable mass transfer might not be generally achievable for binary NSs. 



Law 



32 G 4 MV 



(12) 



5 c 5 a 5 



IV. RESULTS AND DISCUSSION 
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FIG. 6: Left panel: The evolution of a 1.5 M Q SQM star (solid line) is compared with that of a 1.5 M© 
normal neutron star (dashed line) with a 3.5 M Q black hole during merger. The compact stars' masses, 
radii, and outspiral orbital separations during stable mass transfer are shown. (For simplicity, the inspiral 
evolutions are not indicated.) Right panel: Luminosities of gravity wave emission for the two cases. 

Note that moderate mass NSs have a ~ 0, so Eq. requires q < 5/6. Had we used the more 
exact formula of Eggleton, Eq. (ffij). we would have found an even greater restriction q < 0.78. For 
example, the binary pulsar PSR1913+16 has q = 0.96. In contrast, for small to moderate mass 
self-bound SQM stars a ~ 1/3, so that stable mass transfer can occur for q < 1, which is its entire 
domain! 

The evolution of self-bound SQM star - black hole merger including the effects of non- 
conservative mass transfer, tidal synchronization, the presence of an accretion disk, etc., together 
with extensions to include further effects of general relativity, will be reported separately. 

In conclusion, compact star mergers offer a tantalizing possibility for the detection via gravita- 
tional waves of strange quark matter in the form of self-bound stars. 
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